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Abstract 

We have performed a detailed phenomenological investigation of the hard 
colour singlet exchange process which is observed at the Tevatron in events 
which have a large rapidity gap between outgoing jets. We include the effects 
of multiple interactions to obtain a prediction for the gap survival factor. 
Comparing the data on the fraction of gap events with the prediction from 
BFKL pomeron exchange we find agreement provided that a constant value 
of as is used in the BFKL calculation. Moreover, the value of Os is in line 
with that extracted from measurements made at HERA. 



On leave of absence from * 



1 Introduction 



Diffractive scattering has yet to be understood within the framework of QCD. 
The problem arises because diffractive processes usually depend upon physics at 
distances which are not small on the scale typical of strong interactions. 

Diffractive events are often selected by looking for regions in rapidity which are void 
of activity. Rapidity gaps are correlated with diffractive scattering since diffrac- 
tion takes place in the limit where the invariant centre-of-mass energy is much 
larger than any other invariant. The challenge is to understand the gap producing 
mechanism using QCD. 

Figure 1 shows a generic rapidity gap event and serves to set our notation. The 
final state is composed of two clusters, X and Y, which are separated in rapidity. 
As usual the Mandelstam variables are labelled s and t. 




Figure 1: A generic rapidity gap process. 



Typically, the gaps are produced with little or no momentum transfer, t, and this 
makes life difficult for those wanting to employ the tools of perturbative QCD. 
However, it is possible to squeeze the gap producing mechanism to short distances: 
one looks for gap events in which a large momentum is transferred across the gap, 
i.e. -t > A^jcD [1,2]. 

Experiments at the Tevatron and at HERA have already examined a variety of 
these large t diffractive processes, the most inclusive of which is the double dissoci- 
ation process '-/p X Y which has been studied at HERA [3].^ The gaps-between- 
jets process, in which systems X and Y are each required to contain at least one 
jet, have been studied both at the Tevatron [4-7] and at HERA [8,9]. Additionally, 
the vector meson production process, jp ^ V Y , in which the system X is a single 
vector meson, V, has been measured at HERA [10,11]. 

^At HERA system X is the product of the photon dissociation and system Y is the product 
of the proton dissociation. 
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Theoretically, the presence of a large momentum transfer across the gap supports 
the use of perturbative QCD [12]. The non-perturbative physics can be factorised 
into the usual parton density functions which can then be convoluted with the hard 
partonic scattering subprocess to determine the physical cross-section. The hard 
scattering subprocess can be calculated to an accuracy where all leading logarithms 
in energy, i.e. all terms ~ Q;"ln"'(s/|i|), are summed up using the formalism of 
BFKL [13]. Calculations have been performed for all of the above processes to 
leading logarithmic accuracy, see [14] for a recent summary. 

In this paper, we wish to examine the Tevatron data on the gaps-between-jets 
process. We perform a detailed analysis of these data at 630 GeV and 1800 GeV and 
are drawn to conclude that the data are fully consistent with the leading logarithmic 
calculation of BFKL providing one fixes the strong coupling. This conclusion is 
consistent with that which can also be drawn from the HERA data collected to 
date. In our study, we pay particular attention to the issue of gap survival: we use 
Pythia to model the underlying event [15]. We consider observables at both the 
parton and hadron level. 



The theoretical calculation 



We start by explaining how we perform our theoretical calculations. We use the 
hard partonic subprocess first calculated in [2]: 



da{qq qq) 



dt \ TT J (iV2- 1)2^2 

where the leading logarithm BFKL kernel is 



du- 



ly 



(1) 



TT 



Re[^(l) - V;(l/2 + ii/)] 



and 



y = ^V = '^^{^ 



(2) 



(3) 



To this level of accuracy, the jets have equal and opposite transverse momenta, 
p"^ = —t. Throughout this paper we evaluate the right-hand- side of (1) without 
approximation although we do note that, in the limit At; S> 1, it simplifies to 



dt 



{7a,CAmyy 



(4) 



where loq — u;{0) — C^(4 In 2/77)0:5. The gluon-quark and gluon-gluon subprocesses 
are the same as the quark-quark subprocess, up to colour factors, and so we can 
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write the parton level cross-section for the gaps-between-jets process as 



da{hi h2^XY) 
dxidx2dt 



81 



16 

da(qq 



16 



qq) 



dt 



(5) 



where gi{xi, /j,"^) is the gluon parton density function for hadron i and T,i{xi, n"^) 
is the sum over all quark and antiquark density functions for hadron i. These 
sub-processes have been coded into Herwig [16,17] thereby allowing one to com- 
pute the likely corrections due to parton showering and hadronisation. We use a 
particular version of Herwig which includes the BFKL subprocesses obtained us- 
ing (1) without approximation [18].^ Subsequently we take the factorisation scale 
/i^ — —t (in the region of interest, s ^ — t, this is a very good approximation to 
the value chosen in Herwig). We have used both CTEQ2M and CTEQ3M parton 
distribution functions [19,20], and have found the differences to be very small. 
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Figure 2: D0 data compared with a BFKL calculation. Plot from [5]. 



^The code can be obtained on request from coxb@mail.desy.de. 
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D0 data versus the BFKL pomeron 



Let us now summarise the situation to date. Figure 2 shows the D0 data compared 
to a BFKL calculation. What is shown is the fraction of all dijet events in which 
the central region of the detector \rj\ < 1 is devoid of hadronic activity (see later 
for a more detailed explanation of the D0 selection cuts) as a function of the 
transverse energy of the lower Et jet, Et2, and the separation in rapidity of the 
jets, Arj. The CDF data at 1800 GeV do not go down to such low values of Et 
and are consistent with a fiat Et distribution. 

The BFKL curve is quite clearly excluded by the data. However, it should be 
understood that this curve was generated using the default settings in Herwig [17]. 
In particular, this means that the asymptotic cross-section of (4) was used with 
a fixed value of uq = 0.3 in the exponent and a fixed value of as = 0.25 in the 
denominator, whilst the prcfactor was allowed to run with —t according to the 
two-loop beta function. These settings are specific to the BFKL subprocess, in all 
other subprocesses the Herwig default is to use the two-loop running couphng. 
The falling of the BFKL curve with increasing Et2 is driven by the running of the 
coupling in the pre- factor since the gap fraction goes like ~ otg/a^. 

Now, to leading logarithmic accuracy as is simply an unknown parameter. Higher 

order corrections will indeed cause the coupling to run, however it is not clear how 
this should be done in a consistent way. In this paper we restrict ourselves to 
the leading logarithmic approximation and treat the coupling as a free parameter. 
Moreover, we are guided by recent HERA data on the double dissociation process 
[3] which can be described by the leading logarithmic BFKL formalism with a 
fixed strong coupling, as = 0.17. Wc also note that a fixed coupling constant was 
needed in order to explain the high-t data on pp elastic scattering via three gluon 
exchange [21]. Furthermore, NLO corrections suggest a fixed value for the leading 
eigenvalue of the BFKL equation, ci;(0), [22] which in turn suggests the use of a 
fixed coupling in the LLA kernel, i.e. in equation (2). 

Underlying events and gap survival 

So far we have said nothing about the possibility that gaps produced by the colour 
singlet exchange might be filled in as a consequence of secondary interactions be- 
tween the colliding hadrons. The possibility that gaps do not survive formally 
invalidates the factorisation of the cross-section which we assumed in the previous 
sub-section. However, there are fairly good reasons (as we shall see in the next 
section) for expecting the effect to enter as a multiplicative 'gap survival' factor, 
S, which specifies the probability that a gap produced at the hard scatter level 
survives to the hadron level. Note that we do expect S to vary with the centre- 
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of-mass energy of the colliding beams: as the centre-of-mass energy increases the 
corresponding proliferation of low-x partons makes a secondary partonic interac- 
tion more likely. 

In the next section we shall nsc Pythia to extract the gap survival factors relevant 
to the Tevatron operating at 1800 GeV, at 630 GeV and for HERA at a fixed jp 
energy of 200 GeV. We will then use these factors to correct our cross-sections for 
processes involving gaps, i.e. 

dajh h2^XY) ^ ^ dajh h2^XY) 
dxidx2dt dxidx2dt ' 

where S is the fraction of gap events which survive after including the effect of the 
underlying event as modelled in Pythia. 



2 Gap Survival and Multiple Interactions 

To investigate the gap survival probabilities we use a model for multiple interactions 
available in the Pythia program. Here the probability to have several parton- 
parton interactions in the same collision is modelled using perturbative QCD. For 
a given hard parton-parton scattering there is a probability of having additional 
scatterings given by the parton densities and the usual LO 2 — > 2 matrix elements. 
The matrix elements are divergent for p± ^ and are regulated by replacing the 
l/p\ pole with l/ip^ + pIo)^ s-iid using p^ + p5_o argument in rather than 
just p^. The value of p±o varies with collision energy, E^m, and we have used the 
default behaviour: p±o = (^cm/(lTeV))°-i6 • 2.1 GeV. 

The probability for additional interactions is not fixed but varies according to an 
impact-parameter picture, where central collisions are more likely to have multiple 
interactions. The partons in the proton are assumed to be distributed according 
to a double-gaussian as described in [15,23]. There are several parameters in this 
model and we have used the default setting for each.'' This is adequate for the 
accuracy required here but it should be noted that the parameters have not been 
tuned to Tevatron data. 

Hard BFKL pomeron exchange has not been implemented in Pythia yet, but we 
can investigate the effects of multiple interactions on the gap survival probabilities 
in general by looking at high-i photon exchange. In Figure 3 and Figure 4 we show 
the gap fraction in photon exchange events at 1800 GeV with (mi4) and without 
(miO) multiple interactions and, as expected, the probability for a gap introduced 
by the colourless photon scattering is greatly reduced in the presence of multiple 

^Setting the switch MSTP(82)=4 in Pythia, with everything else default, will give the model 
as we have used it. 
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Figure 3: Gap fraction in pho- Figure 4: Gap fraction in pho- 
ton exchange events at 1800 GeV ton exchange events at 1800 GeV 
with/without Mi's: Et2 spectrum with/without MPs: Arj spectrum 



interactions. The reason is that any additional scattering is hkely to be a colourful 
one, destroying the gap. Thus we find the gap survival probability to be a simple 
multiplicative factor almost independent of the gap size and Et although varying 
with the total collision energy. Figure 5 shows the effect of multiple interactions 
at the lower energy of 630 GeV. 

In Figure 6 we show the effect in jp interactions at 200 GeV, as is relevant for 
HERA. In Pythia the procedure for generating these jp events requires some 
care. The cross-section is divided into three parts: direct, anomalous and VMD. 
Direct interactions, where the photon couples directly to the hard interaction, are 
not present in gamma exchange. Anomalous processes are resolved ones where 
the evolution from the photon is purely perturbative, while VMD are resolved 
processes where the photon has fluctuated into a vector meson state. Multiple 
interactions are only included in the latter processes, which is why the gap survival 
probability is much higher than in the hadron-hadron case. For jp we also use a 
simpler approach to multiple interactions, with an impact-parameter-independent 
probability for secondary scatterings since for this option there exists a set of 
tuned parameters [24] which we have used in the generation. The gap fraction 
after including multiple interactions is labelled 'mil' in Figure 6. 

In summary, we find 5(1800 GeV) = 22%, 5(630 GeV) = 35% and 
5(HERA) = 67%. 
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Figure 5: Gap fraction in pho- Figure 6: Gap fraction in photon ex- 
ton exchange events at 630 GeV change events at HERA with/without 
with/without Mi's Mi's 



Note that requiring a large rapidity gap not only selects events with colour- singlet 
exchange in general, but in particular events with colour- singlet exchange without 
additional interactions. 

Multiple interactions also give rise to the so-called jet pedestal and underlying event 
effects. This means that the jets measured in hadron-hadron collisions cannot be 
compared directly to e.g. predictions from fixed order perturbation theory. In 
Figure 7 we show jet profiles obtained from Pythia with and without multiple 
interactions (and with \S(f)\ < 0.7). The proton remnant is at 61] > 0. It is clear 
that multiple interactions introduce a jet pedestal of more than 1 GeV of Et per 
unit rapidity. For comparison, also shown is the jet pedestal from Herwig. We 
note that Herwig predicts a greater amount of energy outside the jet cone than 
Pythia without multiple interactions. On the remnant side, 61] > 0, the difference 
is accounted for by different treatments of the remnant (running Herwig with 
IOPREM=0 reproduces the Pythia result). However, on the gap side, 6ri < 0, the 
large difference is not so easy to explain. Herwig does use a larger coupling than 
Pythia (see later), although this is not enough to account for the differences. The 
disagreement could be related to the differences in treatment of QCD coherence. 
Coherence in Herwig is achieved by ordering emissions in the parton cascades 
in angle, while in Pythia the ordering is in virtuality with an additional angular 
constraint. The double ordering in Pythia reduces the emission probability and 
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Figure 7: Jet rj profiles 



the fact tfiat we are looking at jets at large rapidity (small angles) may enhance this 
effect. We have checked that relaxing the angular constraint in Pythia increases 
the Et-S-Ow on the gap side by ~ 20%, although this needs to be studied further 
before firm conclusions can be made. Fortunately, as we shall see, it is likely that 
these differences mostly influence the overall normalisation of the gap fraction. 
However, it is important that these jet shapes are thoroughly confronted with data 
before precise estimates of the hadronisation corrections can be made. 

In the D0 jet measurements the excess Ex from the underlying event is taken 
into account by correcting the jet E^ using data. In particular, the correction is 
determined by looking at the E^ in regions away from the jets. The correctness of 
this subtraction is also supported by a study of minimum bias events. The bottom 
line is that D0 subtract around 1 GeV from the E^ of each reconstructed jet [25]. 
In particular, in the gap fraction measurement, this subtraction is performed for 
all jets, including those in gap events. But, as we already pointed out, requiring 
a large rapidity gap also selects events without multiple interactions, where the 
jet pedestal is absent, or at least much smaller. Now, 1 GeV may seem like a 
small correction when we are looking at jets of Et = 18 GeV or larger, but recall 
that the Et spectrum is very steeply falling: typically it falls faster than 1/E^. 
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The correction may therefore be as large as 30% for Et — 18 GeV, while for 
Et > 50 GeV it will be less than 10%. Looking at Figure 2 we can roughly undo 
the erroneous correction of the jets in the gap events. The lowest Et bin would 
then be increased a factor (19/18)^ ~ 1.24 while the other points would be less 
affected. The increase of the gap fraction with increasing Et then becomes less 
pronounced. 

3 Results 

Cross-sections at the parton and hadron level 

In this section we show our results for jet cross-sections. Specifically, we work 
at 1800 GeV centre-of-mass energy and make the same analysis cuts as the D0 
collaboration [5]. Jets are found using a cone algorithm [26,27] with cone radius 
0.7 and the OVLIM parameter set to 0.5. The inclusive dijet sample is defined by 
the following cuts: 

• 1^1 1 ) 1^72! > 1-9, i.e. jets are forward or backward 

• ViV2 < 0, i.e. opposite side jets 

• Et2 > 15 GeV 

• A77 > 4, i.e. jets are far apart in rapidity. 

The sub-sample of gap events is obtained by employing the further cut that there 
be no particles emitted in the central region \ri\ < 1 with energy greater than 
300 MeV. The CDF collaboration have performed a similar analysis [7] and we 
comment on their data wherever relevant. GDF use the same selection criteria as 
D0 except that they chose 1.8 < |77i,2| < 3.5, Et2 > 20 GeV for the 1800 GeV 
sample and Et2 > 8 GeV for the 630 GeV sample, and they make no further cut 
on the jet separation.^ 

We note that observables like the fraction of gap events to all events (the so-called 
gap fraction) are strictly not infrared safe. The definition of a gap event requires 
the specification of some threshold energy. Particles emitted into the gap with 
energy below this threshold by definition do not spoil the gap. In theoretical 
calculations this leads to logarithms in the infrared scale which formally diverge as 
the threshold is taken to zero. Thus theoretical predictions based on perturbation 
theory are in principle unstable as the threshold is lowered. One way to avoid this 
problem is to define a 'gap' to be an interval in rapidity which does not contain 

^Strictly speaking we should set 0VLIM=0.75 when comparing to CDF data. However, the 
effect of implementing this change is negligible. 
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any jets above some Ex which is perturbatively large. However, the CDF and D0 
gap samples are defined in such a way as to greatly reduce the sensitivity to the 
threshold energy. By insisting that the gap be a strip of width at least 2 units in 
rapidity centred on = and by forcing the jets to be at > 1.9 the majority of 
events have gaps which do not encroach too closely to the edge of the jets. As a 
result, any soft radiation tends not to be emitted into the gap (it is localised around 
the leading parton). The problem would have been more acute if the gap sample 
had been defined by insisting that no particles be emitted (above the threshold) 
between the edges of the jet cones. This latter protocol has been used by the 
HERA experiments and should lead to results which depend significantly upon the 
threshold energy. We have exphcitly checked that our results for the Tevatron do 
not change significantly as the threshold energy is varied from MeV to 300 MeV. 

In all cases results labelled 'bfkl fixed ctg' represent the cross-sections obtained 
using only the 2^2 hard subprocess cross-section of (1) with = 0.17; the value 
hinted at by the HERA data. Those labelled with 'bfkl running a^' are obtained 
using the running coupling with the two-loop beta function (Aq^d = 200 MeV) 
in the prefactor whilst keeping a fixed coupling, as — 0.17, in the evaluation 
of uj{v). The parton showering and hadronisation are determined by the usual 
Herwig defaults. Note that we have not yet corrected the BFKL data by the gap 
survival factor determined earlier. 

All other results are cross-sections for the non-BFKL 2^2 subprocesses, i.e. they 
correspond to the usual QCD non- colour-singlet exchange subprocesses. More 
specifically: curves labelled 'Herwig lo' are generated by Herwig with a one- 
loop running coupling (Aqcj-, = 362 MeV) whilst those labelled 'Herwig nlo' are 

the Herwig default curves obtained using a two-loop running coupling (A^^^ 
200 MeV); similarly the Pythia curves are generated without multiple interac- 
tions ('miO') or with multiple interactions ('mi4') and correspond to a one-loop 
running coupling with Aqcd determined by the value used in fitting the chosen 
parton density functions (Aq^j-, ^ 200 MeV when using CTEQ2M). The Aqcd 
values quoted are the ones used to generate the parton level theory curves shown 
on some of the plots and, at two loop, are quoted in the MS scheme. The parton 
level theory curves are obtained from those bare 2^2 subprocess cross-sections 
which involve t-channel gluon exchange (t-channel quark exchange gives a negli- 
gible contribution), i.e. they do not include the effects of parton showering and 
hadronisation. 

In Figure 8 we show the Et spectrum of the lower Et jet after hadronisation and 
parton showering. Note that there is a significant difference between the Herwig 
and Pythia predictions. This is due to the different ways of running the coupling 
employed in the two programmes and the differences in fragmentation discussed 
earlier: data would help resolve the issue. Further sensitivity to the choice of as 
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Figure 8: Cross-sections at the hadron level 



(even in the non-BFKL case) can be seen in Figure 9, where we show the Pythia 
Et spectrum of the lower Ej- jet after hadronisation and parton showering but 
without multiple interactions, and compare to the parton level theory predictions 
and to Herwig. Hadronisation and parton showering corrections are clearly very 
significant in both the BFKL and non-BFKL data sets. We have also looked at the 
relative role of hadronisation and parton showering and find that each contributes 
significantly to the deviations from the parton level theory predictions. 

In Figure 10 we compare the parton level theory curves to the parton level cross- 
sections obtained from Monte Carlo. These Arj distributions are obtained for the 
'low Et' jet sample, i.e. for those jets in the range 15 GeV < Et2 < 25 GeV. As 
well as providing a cross-check, they serve to illustrate the dominance of t-channel 
gluon exchange. This plot also illustrates the effect of insisting that |?7i|, |?72| > 1-9. 
The fall of the curve at low Arj is a result of the increasingly strong requirement 
that the incoming partons have nearly equal and opposite energies. As a result, in 
the absence of the At] > 4 cut, these cross-sections would fall smoothly to zero at 
Af] = 3.8. This kinematic effect could be removed by relaxing the constraint that 
the gap be central, i.e. by allowing the whole dijet system to tilt. The fall off of 
the Af] distribution at large Ai] is due to the decrease of the parton densities with 
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Figure 9: Comparison of Monte Carlo hadron level with theory parton level 



increasing x. 



Gap fractions 

We now turn our attention to the gap fraction. This is the quantity usually mea- 
sured by experiment. To construct the gap fraction we form the ratio 



/ 



N'bfkl + N'o 



QCD 



N, 



BFKL 



(7) 



QCD 



where N^pj^^ is the number of events in the BFKL sample which also satisfy the 
gap cut, Nq(jj~, is the number of events in the non-BFKL sample which satisfy the 
gap cut and the sum in the denominator is over all events in the dijet sample. The 
numerator is overwhelmingly dominated by BFKL events in all the plots we show. 

In Figure 11 we show the hadron level gap fraction as a function of Et2- All 
our gap fraction plots include the Pythia gap survival correction factor. The 
stars, corresponding to the ratio 'bfkl fixed — 1 GeV / pythia miO', are the 
result of simulating the D0 jet correction procedure. We have subtracted 1 GeV 
from the BFKL jets and generated the non-BFKL sample by running Pythia 
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Figure 10: Low Et events. Comparison of Monte Carlo parton level with theory 
parton level 

without multiple interactions. Notice that these data points are in agreement with 
the the triangles labelled 'bfkl fixed a^/ pythia mi4' which have been obtained 
without applying the jet correction to the BFKL sample and including multiple 
interaction effects in the non-BFKL sample. This agreement confirms our earlier 
claim that if one should want to correct for the affect of the underlying event (i.e. 
multiple interactions) then the correction can be approximated by a subtraction 
of 1 GeV from each jet. The squares in Figure 11 correspond to the gap fraction 
properly corrected for the underlying event, i.e. only the denominator (non-gap 
sample) has been corrected. Also shown in Figure 11 are comparisons with the 
ratio predicted using Herwig to generate the non-BFKL sample and with a parton 
level theory calculation, labelled 'pythia miO parton' which was computed using 
the same prescription as Pythia for the running coupling (as detailed above). 

The rise at low Et2 is mainly a result of including multiple interactions since, at 
low -E't2, the small amount of energy arising from the underlying event produced 
by secondary scatters enhances significantly the denominator (i.e. the full dijet 
sample) thereby reducing the gap fraction. Although we do note that a weaker 
rise, due to parton showering and hadronisation corrections, is still present in the 



13 



0.05 
0.045 
0.04 



0.035 



0.03 
0.025 
0.02 



0.015 



0.01 
0.005 


15 20 25 30 35 40 45 50 55 

Figure 11: Gap fraction function of Et2 
Herwig and Pythia miO samples. 

Note also the decline of this rise of the gap fraction as Et2 increases. This effect 
is also present in the theory level curve (where it is seen as a falling distribution), 
even for fixed and reflects the tendency for the higher Et2 events to have a 
smaller At] due to the presence of the {xiX2s/p'^)'^'^° factor in the numerator (the 
strong fall off of the parton density functions at large x means that larger pt tends 
to occur without a corresponding increase in mean x). We note that we have used 
the BFKL formalism to sum the leading logarithms in the numerator of the gap 
fraction (supplemented with parton showering and hadronisation as determined 
by Herwig) whilst using established Monte Carlos to predict the denominator. A 
strictly leading logarithmic calculation would also sum the leading logarithms in the 
denominator which could give rise to a factor in the denominator of ~ {xiX2s/p'^)'^° . 
This physics has not yet been incorporated in Herwig or Pythia and would have 
the effect of slowing down the fall of the Et spectrum with increasing Ej'. 

Figure 12 shows the gap fraction as a function of Arj for the low E^ sample and 
Figure 13 shows it for the high E^ sample. Notice that the effect of multiple 
interactions on the denominator of the gap fraction diminishes as Et2 increases. 

Collectively, these gap fraction plots illustrate the uncertainty in theoretical pre- 
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Figure 12: Gap fraction as a function Figure 13: Gap fraction as a function 
of At] for 15 < Et2 < 25 GeV of Ar] for Et2 > 30 GeV 

diction even in the case of a particular ansatz for the BFKL behaviour: parton 
showering and hadronisation corrections are significant in computing the gap frac- 
tion. 



Comparison to data 

Now we come to our comparison with the D0 data. In order to make the com- 
parison, we must first correct our gap events as discussed in Section 2, i.e. we 
subtract 1 GeV from each jet in the BFKL sample and generate the non-BFKL 
sample using Pythia without multiple interactions. 

Figure 14 shows that after fixing the coupling and making the 1 GeV jet correction 
the open stars agree well with the data. Note that we have chosen to renormalise 
our results by a factor 0.6 in those plots where we compare to the D0 data. That 
this is a reasonable thing to do can be appreciated once it is realised that our re- 
sults have not been fitted to the data and that the overall normalisation is acutely 
sensitive to the magnitude of a^. Furthermore, the overall normalisation of the 
BFKL cross-section is uncertain since, within the leading logarithmic approxima- 
tion, we are free to add any small constant to y in (1), i.e. the transformation 
{s/p\Y^ — ^ {s/{cp^)Y^ where c is some unknown constant ~ 1 is perfectly admis- 
sible. Given these points, we conclude that the D0 data are in agreement with 
the leading order BFKL result. We note that it is not possible to obtain agree- 
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Figure 14: Gap fraction compared to the D0 data: Et2 spectrum. 



ment with the data using a strong couphng which runs with the jet E^, even after 
making the 1 GeV subtraction (the sohd circles in Figure 14). 

We could have chosen to use Herwig to model the non-BFKL sample. In which 
case, the only significant effect, as illustrated in the preceding plots, would be a 
further overall shift in normalisation. 

Figure 15 shows the At] distribution at low Et and again the shape agrees well 
with the data. A similar story can be seen in Figure 16 which shows the At] 
distribution at high Et- The data hint at the possibility that a renormalisation 
somewhat closer to unity (compared to the lower Et case) of the BFKL prediction 
is needed in order to bring agreement with the data. However, as we mentioned 
earlier, summing the Inl/x terms in the denominator will reduce the fall of the 
Ej- spectrum as Ej- increases and this would in turn lift the theory predictions in 
Figure 16. 

In Figure 17 we show our prediction for the D0 gap fraction at 630 GeV compared 
to that at 1800 GeV. D0 has measured the ratio of the number of gap events at 
630 GeV to the number at 1800 GeV and finds a value of 3.4 ±1.2 [5].^ This 



^We note that the D0 ratio was actually obtained for Et2 > 12 GeV. 
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Figure 15: Gap fraction compared to Figure 16: Gap fraction compared to 
the D0 data: Ar/ spectrum at low Ex- the D0 data: Arj spectrum at high Ex- 



value can be seen to accord with the results shown in Figure 17. The data points 
represent our prediction whilst the solid line represents the parton level theory 
prediction. Clearly the effects of parton showering and hadronisation are very 
significant. The rise of the ratio illustrated in Figure 17 at large Arj is understood. 
The restriction x < 1 forces the gap and non-gap cross-sections to fall to zero 
at some maximum Ai], Ar/max- Now, the colour connection which exists between 
the jets in the non-gap sample drags the jets closer together in rapidity. This has 
a small effect away from Ar^max (since the Arj spectrum is roughly flat) however 
as At] — >• A?7max it leads to a more rapid vanishing of the non-gap cross-section 
than occurs in the gap cross-section. This effect, combined with the fact that 
Ar7max(630 GeV) < Ar7max(1800 GeV) leads to the rise illustrated in Figure 17. In 
Figure 18 we show explicitly the gap fraction at 630 GeV and one can see the large 
fragmentation corrections which occur as Arj A?7max- We note that CDF has 
also measured the ratio as a function of the momentum fractions of the final state 
jets (which are approximately equal to the momentum fractions of the incoming 
partons). They find a ratio of 2.4 ± 0.9 at fixed momentum fractions (the value is 
independent of the momentum fraction) [7]. The fall of the solid line (representing 
the tree level prediction) in Figure 17 arises since the 630 GeV data are at higher 
parton x than the corresponding 1800 GeV data and are consequently in a region 
where the gluon density falls off more sharply with increasing x than the quark 
distributions. Since the gluon is given more weight in the BFKL sample, this fall 
compensates the rise of the gap fraction due to the exp(2a;o|/) factor leading to a 
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Figure 17: Ratio of the gap fraction at 
630 GeV to that at 1800 GeV: At] spec- 
trum for > 15 GeV. 



Figure 18: Gap fraction at 630 GeV: 
Ar] spectrum for > 15 GeV. 



flatter Arj distribution at 630 GeV than at 1800 GeV. 

Finally, in Figure 19 we compare our results with the CDF data at 630 GeV 
and 1800 GeV: t]* = Ar}/2. Note that CDF do not attempt to correct their jets 
to include the effect of an underlying event. In this plot, our theory points are 
obtained using a renormalisation factor of unity (compared to 0.6 in the D0 case) 
which is more in line with the normalisation suggested by the higher Ex D0 data 
and, as we have already discussed, is possibly expected on theoretical grounds after 
the inclusion of BFKL effects in the denominator of the gap fraction. We then find 
reasonable agreement with the data except at the larger values of rj* where we are 
quite unable to explain a fall in the Vj* distribution. Recall however that D0 do not 
see a fall at large Arj. Further clarification of the situation will require an increase 
in statistics. 



4 Conclusions 

We have explicitly demonstrated that the Tevatron data on the gaps-between-jets 
process at both 630 GeV and 1800 GeV are in broad agreement with the predictions 
obtained using the leading order BFKL formalism. However, we are not able to 
explain the behaviour of the CDF gap fraction at large Ar/. Agreement is obtained 



18 



0.05 

V 

0.045 
0.04 

0.035 
0.03 



0.025 
0.02 

0.015 
0.01 



0.005 





1.8 2 2.2 2.4 2.6 2.8 3 3.2 3.4 

Figure 19: Gap fraction compared to the CDF data at 630 GeV (open squares) 
and 1800 GeV (open circles). 

using the same fixed value of = 0.17 as was used to explain the recent HERA 
data on high-t double diffraction dissociation. 

This agreement should be set in context. The theoretical formalism used does suffer 
from being evaluated only to leading logarithmic accuracy. There is an essentially 
unknown overall normalisation (due to the ambiguity in the choice of scale which 
defines the leading logarithms, i.e. the parameter c discussed in Section 3) and 
the procedure for determining the correct treatment of the running coupling is not 
defined. To improve upon these shortcomings requires an understanding of BFKL 
dynamics at non-zero t which goes beyond the leading logarithmic approximation. 

Crucial in developing a description which is simultaneously consistent with the 
Tevatron data at 630 GeV and at 1800 GeV, and the HERA data at 200 GeV, is 
a realistic treatment of the soft underlying event and its impact on the survival 
of rapidity gaps. Using Pythia to model the underlying event we find an overall 
gap survival factor which is independent of the jet Et and rapidity but dependent 
upon the overall centre-of-mass energy. 

However, the gap fraction is not a very clean observable. We found very significant 
corrections from parton showering and hadronisation to both the colour- singlet 
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and non-colour singlet exchange processes which do not cancel on taking the ratio 
and which persist to the largest values of E^- Moreover, the influence of the soft 
underlying event (which we accounted for using a multiple interactions model in 
Pythia), is also very significant at the smaller values of Et, i.e. Et < 30 GeV. 
A good understanding of the role of the underlying event for Et < 30 GeV is 
therefore crucial in comparing with the data. 

Acknowledgements 

We should hke to thank Andrew Brandt, Dino Gouhanos, Mark Hayes, Mike Sey- 
mour and Torbjorn Sjostrand for helpful discussions. This work was supported 

by the EU Fourth Framework Programme 'Training and Mobility of Researchers', 
Network 'Quantum Chromodynamics and the Deep Structure of Elementary Par- 
ticles', contract FMRX-CT98-0194 (DG 12-MIHT). BC would like to thank the 
UK's Particle Physics and Astronomy Research Council for support. 

References 

[1] J. D. Bjorken, Phys. Rev. D47 (1992) 101. 

[2] A. H. Mueller and W. -K. Tang, Phys. Lett. B284 (1992) 123. 

[3] HI Collaboration, "Double Diffraction Dissociation at large \t\ in Photopro- 
duction at HERA" , contribution to the 29th International Conference on High- 
Energy Physics ICHEP98, Vancouver, Canada, 1998; 

B. E. Cox, "Double Diffraction Dissociation at large \t\ from HI", contribution 
to the DIS99 Workshop, Zeuthen, Germany (1999) hep-ph/9906203. 

[4] S. Abachi et al (D0 Collaboration), Phys. Rev. Lett. 72 (1994) 2332; 
Phys. Rev. Lett. 76 (1996) 734. 

[5] B. Abbott et al (D0 Collaboration), Phys. Lett. B81 (1998) 189. 

[6] F. Abe at al (CDF Collaboration), Phys. Rev. Lett. 74 (1995) 855; 
Phys. Rev. Lett. 80 (1998) 1156. 

[7] F. Abe at al (CDF Collaboration), Phys. Rev. Lett. 81 (1998) 5278. 

[8] M. Derrick et al (ZEUS Collaboration), Phys. Lett. B369 (1996) 55. 

[9] HI Collaboration, "Rapidity gaps between jets in Photoproduction at HERA" , 
contribution to the International Europhysics Conference on High Energy 
Physics, Jerusalem, Israel (1997). 
[10] HI Collaboration, "Production of J/ip Mesons with large \t\ at HERA" , contri- 
bution to the International Europhysics Conference on High Energy Physics, 
Jerusalem, Israel (1997). 



20 



[11] J. Crittenden, "Recent Results from Decay-Angle Analyses of p° Photopro- 

duction at High Momentum Transfer from ZEUS" , contribution to the DIS99 

Workshop, Zeuthen, Germany (1999) hep-ex/9906005. 
[12] J.R. Forshaw and RJ. Sutton, Euro. Phys. J. C14 (1998) 285. 
[13] I. Balitsky and L.N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822. 
[14] J. R. Forshaw, "High-t Diffraction", contribution to the DIS99 Workshop, 

Zeuthen, Germany (1999) GERN-TH/99-143, hep-ph/9905557. 
[15] Pythia version 6.127, program and manual, 

T. Sjostrand, Comput. Phys. Gomm. 82 (1994) 74. 

See also http://www.thep.lu.se/~torbjorn/Pythia.html. 
[16] G. Marchesini et al, Gomp. Phys. Gomm. 67 (1992) 465. 
[17] J. M. Butterworth, M. E. Hayes, M. H. Seymour and L. E. Sinclair, "Rapidity 

gaps between jets" , in the proceedings of the Workshop 'Future Physics at 

HERA', eds. G. Ingelman, A. de Roeck and R. Klanner, DESY (1996). 

[18] M. E. Hayes, University of Bristol PhD Thesis (1997). 
[19] H. L. Lai et al, Phys. Rev. D55 (1997) 1280. 

[20] H. Plothow-Besch, "PDFLIB User's Manual", W5051 PDFLIB, 1997.07.02, 
GERN-PPE; Int. J. Mod. Phys. AID (1995) 2901. 

[21] A. Donnachie and P. V. Landshoff, Z. Phys. C2 (1979) 55, erratum-ibid C2 
(1979) 372; Phys. Lett. B387 (1996) 637. 

[22] S. J. Brodsky et al, "The QGD Pomeron with Optimal Renormalisation" , 
SLAG-PUB-8037, IITAP-98-010, hep-ph/9901229. 

[23] T. Sjostrand and M. van Zijl, Phys. Rev. D36 (1987) 2019 

[24] J.M. Butterworth and R.J. Taylor, "A Global study of Photon In- 
duced Jet Production", hep-ph/9907394, to be published in the 
proceedings of Photon '99, Freiburg, Germany (1999). See also 
http : / / www . hep . ucl . ac . uk/~ jmb/HZTOOL/. 

[25] B. Abbott et al (D0 GoUaboration) , Nucl. lustrum. Methods A424 (1999) 
352. See also http://www-dO.fnal.gov/~jkraiie/approved_plots. 

[26] S.D. Ellis, private communication to the OPAL GoUaboration; D.E. Soper 
and H.-G. Yang, private communication to the OPAL GoUaboration; L.A. del 
Pozo, University of Cambridge PhD thesis, RALT-002 (1993); R. Akers et al 
(OPAL Collaboration), Z. Phys. C63 (1994) 197. 

[27] F. Abe et al (CDF Collaboration), Phys. Rev. D45 (1992) 1448. 



21 



